Distribution of spectral widths and preponderance of spin-0 ground states in nuclei 
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We use a single j-shell model with random two-body interactions to derive closed expressions 
for the distribution of and the correlations between spectral widths of different spins. This task is 
facilitated by introducing two-body operators whose squared spectral widths sum up to the squared 
spectral width of the random Hamiltonian. The spin-0 width is characterized by a relatively large 
average value and small fluctuations while the width of maximum spin has the largest average and 
the largest fluctuations. The approximate proportionality between widths and spectral radii explains 
the preponderance of spin-0 ground states. 
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Introduction: Many regular features in the low-lying 
parts of nuclear spectra can be attributed to a short- 
range, attractive effective interaction with pairing and 
quadruple components. It thus came as a surprise when 
Johnson, Bertsch and Dean d] found that for even-even 
nuclei, an ensemble of nuclear shell-model Hamiltonians 
with random two-body interactions is likely to yield a 
spin-zero ground state. This is especially so since the 
probability for a spin-zero ground state was found to be 
much larger than the fraction of spin-zero states in the 
model space. Subsequent work showed that similar regu- 
larities exist in bosonic [2j and electronic [2j many-body 
systems with random two-body interactions. Thus the 
phenomenon of spin-0 preponderance seems a very ro- 
bust and rather generic feature. 

The phenomenon has received intense attention, with 
reviews in Refs. 0, 0]. Quantitative explanations for 
the spin-0 ground state dominance were given for ex- 
actly solvable boson and fermion 0,0 systems, while 
mean-field theory provides an understanding for the in- 
teracting boson model 0. For more complex fermion 
systems, the situation seems more difficult. Approaches 
based on the ensemble-av erag ed s pec tral widths, while 
useful for certain systems 0, H2I fl3l Ibl , failed to 
provide an explanation for a single j-shell 0, 0, E3 • 
Recently, Mulhall et al. found that the spectral cen- 
troids form a band or an inverted band, depending on 
the sign of the effective moment of inertia. However, it 
is not clear how to relate this result to the spin of the 
ground state. In another approach, Zhao et al. ^3] made 
quantitative predictions based on the diagonalization of 
the individual two-body operators. However, the success 
of this approach is not well understood. 

As a key to the understanding of the preponderance of 
spin-0 ground states, we propose in this Letter the dis- 
tribution of and the correlations between the widths of 
the densities of states with spin J (in short: the spec- 
tral widths). As a model system, we consider n = 6 
fcrmions in a shell with spin j = 19/2 which interact via 



a random two-body interaction. The data points in the 
top part of Fig. ^show the probability that the ground 
state of this system has spin J. This probability is rather 
large for minimum and maximum spin while the relative 
dimensions of the corresponding Hilbert spaces are very 
small. Thus, the figure exhibits the puzzle pointed out in 
Ref. 0]. The solid line in the upper part of Fig. ^ shows 
the probability that spin J has the largest spectral width. 
This probability is sizable only for minimum and max- 
imum spin. Thus, the preponderance of spin-0 ground 
states is suddenly not a surprise any more! This obser- 
vation prompts us to study the distribution functions for 
the spectral widths, and the correlations between spectral 
widths of different spins. This will eventually lead us to 
a semi-quantitative understanding of the preponderance 
of spin-0 ground states. 

The two-body random ensemble (TBRE). We consider 
the Hamiltonian matrix H(J) in the <ij-dimensional 
Hilbert space Ti-(J) of n-fermion states with spin J for a 
given two-body interaction. The latter has reduced ma- 
trix elements (TBME) v a ,a = 1, ...,a. In our model 
system, there are a — 10 TBME since two fermions can 
have spins j a — 0, 2,4, . . . , 2j — 1. The matrix H(J) is 
linear in the TBME, 

a 

H(J) = 5> Q C Q (J) . (1) 

a 

The matrices C a (J) transport the two-body interaction 
into the Hilbert space TC(J). These matrices are de- 
termined entirely by the geometry of the shell model. 
They are built of coefficients of fractional parentage and 
angular-momentum coupling coefficients. 

In order to obtain generic results, we assume that 
the TBME are uncorrelated Gaussian-distributed ran- 
dom variables with zero mean and unit variance, = 
and v a vp — 8 a p where the bar denotes the ensemble av- 
erage. Our results then apply to almost all two-body 
interactions, the integration measure being the product 
of the differentials of all the v a 's. The matrices H(J) 
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FIG. 1: (Color online) 6 fermions in a shell with spin j — 
19/2. Top: Probability that the ground state has spin J (data 
points) ; probability that spin J has the largest spectral width 
(solid line); probability that the product r.jaj is maximal 
(dashed line). Bottom: Scaling factor rj between the widths 
and spectral radii. Inset: Spectral radius _Ro versus width ao 
(data points) and the linear fit (line) for total spin J = 0. 
(Results from 900 random realizations). 

are sums of random variables and, thus, form a random- 
matrix ensemble, the two-body random ensemble. For 
spins J ^ J', the matrices H{J) and H(J') depend on 
the same random variables v a and are thus correlated. 

Distribution of spectral widths. The spectral width a,j 
of the states with spin J, a random variable, is defined 
as 

a 2 j = dj 1 Tv{H 2 (J)) =J2 v ^ a /3Vp . (2) 

a/3 

Here S a p(J) = dj 1 Ti(C a (J)C/3(J)) is the symmetric and 
positive semi-definite a-dimensional overlap matrix in the 
space of the random variables v a . Diagonalization of 
S = Us 2 U T yields the eigenvalues s\ (J) > s|(J) > . . . > 
s 2 (J) > and the orthogonal, J-dependent matrix U. 
The matrices-in-Hilbert-space B a (J) = J2pUi3 a Cf3(J) 
are orthogonal to each other in the sense of the trace, 
i.e., d~j lr TY{B a {J)Bp{J)) = 5 a0 s 2 a {J), and B a has the 
spectral width s a . In terms of the _Bs, the Hamiltonian 
reads 

H{J) = Y J W a B a {J) . (3) 

OL 

The new Gaussian random variables w a = ^pvpUpa 
obey = and uvuJg = 5 a p . For fixed J, the random- 
matrix model (0) is equivalent to the random-matrix 
model (01. 

In the new basis, we have 

- 2 / = E^ s '( J )- ( 4 ) 

a 
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FIG. 2: (Color online) Roots of the ten eigenvalues of the 
eigenvectors B a (J) versus J for 6 fermions in a j = 19/2 
shell. 



The geometric aspects of the shell model are contained in 
the roots s a > of the eigenvalues. Figure |21 shows how 
the s a vary with spin J for our model system. For each 
spin J there is one particularly large root s\ . This root in- 
creases almost monotonically with J. For each spin there 
is one basis vector B a (J) with root zero. This vector is 
given by the matrix representation of J 2 — J( J+ 1) where 
J denotes the spin operator. Indeed, J 2 — J(J+ 1) anni- 
hilates states with spin J and is a scalar two-body oper- 
ator. Hence, the matrix representation of J 2 — J(J + 1) 
is a linear combination of the matrices C a (J) and, at the 
same time, an eigenvector of S a p with eigenvalue zero. 
For large values of J, there is more than one zero eigen- 
value. This is because here the number of independent 
matrix elements is smaller than the number of TBME. 

We first use Fig. [21 for a qualitative discussion of the 
distribution of widths in the TBRE. Clearly, a 2 = J2 a s a 
attains its maximum value for maximum spin J — J ma x- 
Among the low spins, a 2 dominates because the root si 
is almost constant for low spins, and the remaining non- 
zero roots are exceptionally large for J = 0. As for the 
fluctuations, we use two limiting cases for orientation: (i) 
All roots are equal to s 2 . Then a 2 — as 2 and the r.m.s. 
variance is \/2as 2 . (ii) Only the root s\ differs from zero. 
Then a 2 — s\ and the r.m.s. variance is V2sl Thus, we 
expect that in relation to the average width, the r.m.s. 
fluctuations of ao are smaller than those of the o~j for 
large J. 

For a quantitative analysis, we define the probability 
distribution Pj(cr) = 2a8(a 2 — a 2 ) for finding a value of 
the width between a and a + da. The factor 2a stems 
from the differential da 2 = 2ada. For J = J max there is 
only one non-zero root si, and Pj max is a Gaussian with 
width si. For general J, all the integrations over the w a 's 
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FIG. 3: (Color online) Normalized probability distribution 
functions Pj(cr) versus spectral width a. The histograms re- 
sult from 900 random realizations for 6 fermions in a j — 19/2 
shell. The line shows the theoretical result 0. 



can be done, and the expression for Pj(cr) reduces to 



PA*) = - 

7T 



■ -k arctan 2ts„ 



Ate 1 



1/4 



(5) 



Figure 13 shows Pj(u) for several spins J. This agrees 
very well with the histograms from 900 realizations of 
our TBRE. 

Correlations of widths. The correlations between the 
Hamiltonian matrices H(J) and H(J') induce correla- 
tions between the widths oj and oj* . This fact is borne 
out by the numerical simulation: Under the assumption 
that width correlations can be neglected, we can use the 
distributions Pj(cr) to compute the probability that spin 
J has a larger width than spin J'. Such a calculation 
yields, e.g., that oo > 04 with 60% probability. The 
numerical simulation of the TBRE shows, however, that 
do > 04 for 93% of all realizations. Hence, there must be 
strong correlations between the widths. 

We again begin with a qualitative analysis. Let \uj) 
denote that eigenvector of the overlap matrix S a p(J) 
which corresponds to the largest root si(J). Figure 01 
shows the scalar products (uj\uji) as functions of J — J'. 
Clearly, the eigenvectors \uj) depend only very weakly 
on J and are almost identical for total spins that differ 
by just a few units. In particular, (uj\uji) decreases very 
slowly with J for J' = (upper line in Fig. while the 
decrease is faster for J = J max (lower line in Fig. 0}. 
We see that for low spins up to values of 10 to 20, the 
eigenvector of the largest root and the largest root itself 
are almost independent of J. Because of these strong 
correlations among roots of low spins, the contributions 
of the largest roots si to the widths are simultaneously 
small or simultaneously large for most realizations of the 
TBRE. In this way the correlations significantly enhance 




FIG. 4: (Color online) Data points: Values of the scalar prod- 
ucts (uj\uj/) of the largest-width eigenvectors \uj) of the 
overlap matrices S a0 = dJ 1 Tr(C a ( J)C^ ( J)) for J > ,/'. The 
upper and lower line connect the data points corresponding 
to J' = and J = J max , respectively. 



the probability that spin J = has the largest total 
width among all low spins. Similar arguments show that 
the maximum spin J = J max is most likely to exhibit 
the largest width for high spins. These statements are in 
keeping with the numerical results presented in the upper 
part of Fig. ^ 

For quantitative results, we define the probability 
p( J, J') = (a j — a j,) that spin J has a larger width 
than spin J' , where O denotes the unit step function. 
Introducing the eigenvalues q a of the matrix S a p{J) — 
S a f3(J')i w e can perform the ensemble average. Employ- 
ing integral transforms, we arrive at 



p(J, J') 



dt 



sin J2p arctan 2tqp 



1/4 



(6) 



This expression accurately describes the results of our 
numerical simulation. 

From widths to ground states. The bridge between 
the spectral widths discussed so far and the spin of the 
ground state is spanned by the scaling factors r j . These 
are defined as follows. For each value of J, every realiza- 
tion of the TBRE yields a spectrum Ej^, i — 1, . . . , dj. 
Let Rj — maxj \Eji\ be the spectral radius for spin J. 
For each realization, Rj is given with equal probability 
either by the energy of the lowest-lying level or by that of 
the highest-lying level with spin J. The spectral radius 
is related to the spectral width aj of that realization by 
the inequality Rj > aj. However, our numerical results 
show that the much stronger relation 



Rj ~ rj aj 



(7) 



is approximately valid, with r j constant (independent of 
the realization considered). This is shown in the inset of 
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Fig. ^ The scaling factor r j > 1 (shown in the lower 
part of Fig. encodes information contained in the tails 
of the spectral density. The average level density typ- 
ically decreases exponentially fast for energies close to 
the spectral radius. Thus, rj is expected to exhibit only 
a logarithmically weak dependence on the dimension dj 
and on the average width. We neglect this dependence as 
well as the difference between the absolute values of the 
energies of the highest and the lowest level in the spec- 
trum of each realization. Thus, we identify the product 
Tjcrj with the energy of the lowest-lying state with spin 
J. It is then tempting to assume that the probability dis- 
tribution of Rj is simply given by that of oj scaled by the 
factor rj. Our numerical results show, however, that the 
two distributions, although similar, do not really coin- 
cide. However, the relation Q is sufficiently accurate to 
determine reasonably reliably the probability that rjoj 
is maximal. Indeed, this probability is plotted as the 
dashed line in the upper part of Fig. ^ The agreement 
with the data points is very satisfactory. 

We have obtained similar results for n = 8 fermions 
in a j = 19/2 shell. For the odd-number system with 
n = 7, we also found that minimum and maximum spins 
are most likely to exhibit the largest widths. However, 
other spins also have smaller but sizable probabilities of 
having the largest widths. 

Summary. We have studied the spin of the ground 
state for the shell model with random two-body interac- 
tions. To this end, we have derived closed expressions 
for both, the width distribution functions and the width 
correlation functions. Using these results, we have shown 
that spin-0 and maximum spin are most likely to exhibit 
the largest widths. The spin-0 width is characterized 
by a relatively large average value and rather small fluc- 
tuations, while the maximum spin displays the largest 
average and the largest fluctuations. We have numeri- 
cally established an approximate proportionality between 
spectral widths and spectral radii. This relation is suf- 
ficiently reliable to explain the preponderance of spin-0 
ground states in the shell model. 

This work has been confined to a single shell with spin 
j = 19/2. However, our results depend only upon geo- 
metrical properties of that model and are thus expected 
to apply similarly in other shells. We speculate that sim- 
ilar considerations would also apply to other many-body 
systems. 

The distribution functions were calculated using super- 
positions of two-body operators whose squared widths 
sum up to the average spectral width of the random 
Hamiltonian. These operators may have further appli- 
cations in nuclear structure calculations: Effective inter- 
actions are often obtained from fitting two-body matrix 
elements [l£( . These procedures ideally should avoid gen- 



erating linear combinations of two-body operators with 
small or zero spectral widths. 

Our work poses the intriguing analytical problem to 
derive the eigenvalues s a (J). These are determined en- 
tirely by the geometry of the shell model. Hence, the val- 
ues of the s a (J) should be accessible by group-theoretical 
techniques. An understanding of the scaling factors rj 
seems less difficult. It would require the knowledge of the 
shapes of the average spectra. This would allow us to de- 
termine r j from the position of the lowest level (found by 
integrating the normalized spectrum from — oo to where 
the integral equals 1/2) and the width of that spectrum. 
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